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Abstract 

The structure of the /3-function of massless QCD is considered in an¬ 
alytic perturbation theory (APT). The two-loop analytic effective cou¬ 
pling constant is written in terms of the Lambert W function. The 
method of APT is applied to the gluon propagator in the Landau gauge. 

It is shown that there is an nonperturbative ambiguity in determination 
of the anomalous dimension function of the gluon field. One possible 
resolution of the ambiguity is presented. 


1.Introduction 


Consistent theoretical determination of the effective conpling constant has become 
one of the most important tasks of QCD. A serions difhcnlty is the occurrence of Landau 
singularities in the effective coupling constant in perturbation theory (PT). It is widely 
believed that in realistic quantum held theories these unphysical singularities are artifact 
of PT In QCD they appear at small space-like momenta. For this reason PT 

becomes useless for low-energy calculations. Supposedly, these infrared singularities signal 
that the perturbative vacuum differs from the true vacuum of the theory |0. In recent 
years more and more attention has been focused on this problem in connection with large 
order behavior of QCD perturbative expansions. 

On the other hand, most of the suggestions of the conhnement mechanism include 
assumptions about the infrared behaviour of the effective coupling constant. Thus in the 
popular linear conhnement picture singular infrared behaviour is assumed for the effective 
coupling constant as{q^) ~ q~‘^- It was conhrmed that this behaviour is consistent with 
the Schwinger-Dyson equations (SDEs) of QCD p, |[. 

Independently, in the literature the hypothesis about infrared hnite effective coupling 
constant p has become popular. Several convincing arguments have been proposed to 
support the idea that the running coupling ’’freezes” at low energies p. The conhnement 
mechanism which may be consistent with this hypothesis has been proposed in Refs. 0. 
A linearly rising potential is not assumed in this framework. 

One possible resolution of the ’’Landau-ghost” problem has been suggested long time 
ago in Refs. Bi. where the ’’analytization procedure” was elaborated in PT. Recently, 
this method has been successfully applied to the QCD running coupling @, I, 0, 0, 0, 
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131, 113 . An Important observation has been made in Refs. It was found that the 


modihed solution for as{—q'^) has a hnite infrared limit ds(0), which is stable with respect 
to higher order corrections. The method has been applied to the inclusive decay of a r 
lepton [^], and it was shown that the renormalization scheme ambiguity of the results 
is considerably reduced in APT [^, Furthermore, APT provides a self-consistent 
determination of the running coupling constant in the time-like region 
similar in spirit has been followed by the authors of Refs. 


An approach 


15, 16 


The organization of the paper is as follows. In Sec.2 we present the notation and discuss 
some technical aspects of analytic perturbation theory (APT). The structure of the (3- 
function is studied in APT. We write the two-loop analytic effective coupling constant in 
terms of Lambert W function [|^. Thus we avoid Newton’s iteration method of solution. 
The iterative solution is compared numerically with the exact solution in the infrared 
domain. In Sec.3 we apply the method to the gluon propagator in the Landau gauge. We 
have found that the Kallen-Lehmann analyticity and the RG alone are not sufficient to 
uniquely determine the nonperturbative part of the anomalous dimension function for the 
gluon held. Two different dehnitions for the anomalous dimension function are considered. 
We discuss a possible resolution of the ambiguity. Sec. 4 contains the concluding remarks. 


2. The ’’analytization procedure” for the QCD running coupling 

The /3-function in QCD is given by 

= ( 1 ) 

where is the renormalized coupling constant, and qt is the renormalization 

point. For simplicity, we have assumed that the quarks are massless. We have in mind 
some renormalization scheme, which we do not need to specify at this stage. For con¬ 
venience we choose the covariant Landau gauge. Eq. o dehnes the /3-function order by 
order in PT. However, one can assume Eq. (|^ beyond of PT. The assumption is that the 
non trivial /3 function exists with perturbative asymptotic expansion in powers of a* 

f3{as) ^- — {as) ( 2 ) 

The hrst two coefficients of the series are independent of the chosen renormalization 
scheme. Their values are given by 

/3o = (11 - 2Nf/3), /3i = (102 - 38iV;/3), (3) 

with Nf being the number of flavors. The effective coupling constant is dehned by 

k‘^:^as{-k‘^) = P{as{-k^)), (4) 
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as(yU^) = as- (5) 

Convenient definition of the QCD scale parameter A has been given in Ref.[]^. The 

relation between conventionally defined parameter and A is given by A^n = 

with T] = Pi/Pq. With this dehnition and with normalization (Q) the solution of Eq. is 


In 



47r Pi f 47r/3o \ 

PqUs /3o V PiAs) 


+ '0(«s), 


( 6 ) 


where 

i>{as) = f {1/P{x) - 1/P^^\x))dx, 

Jo 

and P^'^px) denotes the two-loop /3-function. To express as{—kp explicitly as a function 
of beyond the one-loop order, one must solve the transcendental equation (P). Usually, 
iteration procedure is used. 

Let a^jP\—kp be nth-order perturbative solution for ag determined by Eq. (||). Then 
the corresponding analytic solution is given by the Kallen-Lehmann integral 


a 


tP-k^) = - 


p^'^pa) 


TT Jo {a — k"^ — iO) 


da, 


(7) 


where spectral function p^'^pa) is the discontinuity of a^p—kp along the positive k‘^ axis. 
The main result of Refs.[§] is that the value Q;(an)(0) in APT is determined by one-loop 
calculation: 

4;’(0) = fiS(O) = 4ir/ft = 47r/(ll - 2/V//3). 

Using Cauchy’s theorem one can represent the function in the equivalent form 

= + ( 8 ) 

the ’’nonperturbative” term compensates the unphysical contributions in We 

shall assume that 9^^'^ satishes the representation 


e^'^p-kp 


1 , 

TT Jkl {a — kJ — iO) ’ 


(9) 


where k\ = —C(„)A^, and C(„) is the positive real number calculable in PT, for example 
C(i) = C( 2 ) = 1. To third order, in the MS scheme and for Nj = 3, we hnd that C( 3 ) = 1.87 
0. The normalization condition (j^) now reads 


OLn 




( 10 ) 


^Note that the one-iteration solution of Eq.(H) develops unphysical complex singularities for n > 3. 
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Let us consider the structure of the /3-function in APT. From Eqs. (^) and ([I0|) we have 






(a -h 


p^'^\a)da 



e(s + I/O' 


ds, 


( 11 ) 


where and ^ denotes the dimensionless variable ^ = exp{—(j)) = Is?/p^. 

Note that ~ exp(—47r//3oas) for ^ 0. From Eq.([IID we see that /3(as) = /3(.^) is 
analytical in the cut-.^ plane. Alternatively, from Eq.(^) we hnd 



das/dlnp^ = + d‘^'^\p^) 


+ ( 12 ) 


where (d^\x) = — YJk=l / (dvr)^^^, and We can rewrite Eq(p/2[) 

in the form 

(danip^s) (dptips) T (dnpips^ d^ O') ^ (13) 

here, for convenience we have used condensed notations suppressing the superscript (n) 
and arguments to the expressions above {6 = ets.). The nonperturbative piece 

(dnp can be written as a power series for 


n 

(dnp{as.e,e) = Y.Bk{eM(^s)\ 


k=0 


(14) 


where 


B, = e, = B,{e) = 


(-i)"a%(-0) , [dk- 


k\ 


dO^ 


+ 


(47r; 


k+l 


for k>2, (15) 


From the representation (^ we see that 9{p^) and 0{p^) are analytical functions in the 
cut ^-plane. The cut is along the positive real ^-axis ^ > l/c(„). Then, from (|l^ it 
follows that the coefficients Bk{0) = Bk{f,) are analytical in the domain f < l/c(„). In the 
weak-coupling limit they vanish exponentially BkiO ~ ^ ~ ^ detailed study 

of the one-loop /3-function can be found in Ref. [§]. 

Now, let us consider the two-loop running coupling in APT more detailed. It is 
convenient to introduce quantities 


4:71 

To second order Eq. (||) reads 


(—fc^) = ^as{—k‘^) and fd{as) = 


da. 


d\n{p‘^ 


(16) 


ln(2:) = 


viz) 


7] In 1 -|- 


1 1 


py{z) 


(17) 


where we denote = —/c^/A^, y{z) = af\—k‘^) and p = fdi/fdo^ {p = 64/81 for Nj = 
3). We shall accept the flavor condition Nf < 8.05, i.e. /3i > 0. In this case the 
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Landau singularity is occurred in the two-loop running coupling [^, Q . The transcendental 
equation is solved by the iteration method. The solution given by one iteration is 


yit{z) = 


In(^) -h r^ln (l -h ^ ln(^)) 


(18) 


The function ( p!8D provides good approximation, for a small enough value of as- Consider 
(^) as an analytic function in the whole complex z-plane. It has the Landau pole at 
z = 1, the logarithmic branch point sX z = exp{—r]) and standard branch point at z = 0. 
The corresponding analytic solution is given by the Kallen-Lehmann integral (see Refs.^l) 


Vitj 


Az) = -[ 

TT Jo 


= i r ^ 

0 X + z TT J-oo + z 




(19) 


where p^\t) = p^\e^) = Pu\e^^^) = Li/{Ll + Ll) with 

Li = TT + pdj:ccos{{p+ t)/r), L 2 = t + p\xi{r/p), r = \j{p + t)^ + vr^. 

Using the Lagrange Inversion Formula we see from (^) that z=l is not a pole but the 


second order brunch point: y{z) 


for ^ 


Thus, (P!8[) violates the analytical 


properties of y{z) near the point z=l. Therefore, it is not unreasonable to solve exactly 
Eq. ([T7|) and to estimate accuracy of the solution (^). Note that the transcendental 
Eq. (P^D is exactly solvable. Indeed, it is easy to convince that the solution is 


= -z 


pi + w{cy 


where 


c = 


( 20 ) 


( 21 ) 


and W denotes the multivalued inverse of C = hFexp(hF). That is the Lambert W 
function. A detailed review of the properties and applications of this interesting function 
can be found in |^. In this article we use dehnition and notation of branches of W 
according the computer algebra system Maple Release 2 0. Now our task is to choose 
the suitable branches of W in Eq. (p0|). The solution must be regular in the cut complex 
z-plane. Consider the function (|2l|) for 0 < r/ < 1. It maps the cut z-plane onto the three- 


sheet Riemannian surface of C- Combining above consideration with known properties of 
the branches of W we hnd 


hF(0,C) if \arg{z)\ < pTi 
W{C) = <( hF(—1,C) if < o.rg{z) < tt 

1U(1,(C) if —TT < arg{z) < —yx 


( 22 ) 


^In Ref. [13] different conventions are accepted. 
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Taking into account Eqs. (|^, ( pOD and (p^), we write the Kallen-Lehmann representation 
for the corresponding causal solution 



where and 

(iTHFibr))) ’ ‘;t<)=exp(-i-l+i(i-l).). (24) 

Equivalently, using the Cauchy formula we can rewrite Eq. (^), for 0 < 2 ; < 1, as follows 


Va 


iz) = —Re 
V 


1 + 1E(0, 


ez ^ 


rjTT 


-PV 


-Im 


lo X — z 


i + iy(o, 


dx, 


ex ^ 


where PV denotes the principal value of the integral along the real axis from 0 to 1. 
Numerical results for the exact two-loop function (PB| ) as well as for the iterative solution 
( 0 ) are summarized in the Table 1. The relative error for the solution ( 0 ), is about 3% 


-k^/N 




-k^jN 



0 

1.000 

1.000 


6 

0.252 

0.261 

0.50 

0.397 

0.405 


7 

0.245 

0.253 

1 

0.352 

0.362 


8 

0.238 

0.247 

2 

0.310 

0.320 


9 

0.233 

0.241 

3 

0.287 

0.297 


10 

0.229 

0.237 

4 

0.272 

0.281 


15 

0.212 

0.219 

5 

0.261 

0.270 


20 

0.201 

0.208 


Table 1: One-iteration approximation vs exact solution for Nf = 3 

for the considered interval. This gives 8.3% error in the value of the QCD scale parameter 
for Nf = 3 (see Table 2). Following Refs.|^ we may use the average 

1 

MQ) = n ^s{p^)dp 

Q Jo 

for comparison. For Q = 2 GeV we hnd that iterative solution ([^) gives answer to much 
better than 0.5% accuracy. The reference values of aaniMr) are taken from Refs.|§]: 
aan{MT-) = 0.36 ± 0.02 with = 1.78 GeV. It should be noted that in the recent 


article [jT^ the transcendental Eq.([l7D for the running coupling constant has been solved 
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aan{Mr) 

0.34 

0.36 

0.38 

Al?’(MeV) 

610 

710 

820 

A(2)(Mev) 

670 

770 

890 

A^\2GeV) 

0.477 

0.500 

0.523 

A{2GeV) 

0.480 

0.502 

0.525 


Table 2: 


numerically on the physical cut. Nevertheless, exact expression (p3D is itself valuable if 
we have in mind to explore global analytical properties of the solution. 

The procedure is not unique. Indeed, one can redehne the analytic charge according 


) = aan{-k ) +6aani-k ) with 6aan{-k ) 




Uj{cr) 


a 


k^ — fO 


da 


where dcian is a ’’genuine” nonperturbative correction, which cannot be reproduced in 
APT, a;(cr) is a integrable function and K{as,fi) is invariant under the RG. New analytic 
charge is asymptotically free and in the weak limit it reproduces correct perturbative 
series. In principle one can £x this ambiguity using appropriate truncated set of the SDEs. 
Indeed, in the recent publications |T^ reasons have given to believe that Saan cannot be 


zero, and an infrared regular QCD running coupling constant cannot be consistent with 
the SDEs for the gluon propagator. Ghost-free axial gauge has been chosen in this work. 
On the contrary, in Ref. the infrared hnite effective coupling constant has been realized 
within the framework of the SDEs in the Landau gauge. However, it was shown that the 
solution for the gluon propagator does not satisfy the Kallen-Lehmann analyticity. We 
see that situation is not so clear. 

Nevertheless, APT as it was formulated in Refs. ii has advantages. It is a minimal 
causal extension of PT and there are no extra parameters in addition to ag. Another 
attractive feature is the stability of the results in APT [Q. APT may provide a valuable 
approximation in the infrared domain where conventional PT fails. 


3. The gluon propagator in APT 


The method discussed above originally has been presented for calculating modihed 
propagators within framework of PT [^. However, the author of |]l[ did not explore 
renormalization invariance of the theory. This step was taken in paper |^, where the 
procedure has been formulated for the effective coupling constants in renormalizable the¬ 
ories. The method was used for the photon propagator. In QED due to the Ward identity 
simple relation holds 

= P{a)/a, (25) 

here 7 ^ (a) is the anomalous dimension of the gauge held. Because of (pS]) the QED 
effective coupling is proportional to the transverse dimensionless structure of the photon 
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propagator. In QCD the relation (p5[) generally is not valid ^ Therefore, the method 
can not be applied to the gluon propagator unambiguously . Let us consider the gluon 
propagator in the Landau gauge D^y{k) = {g^j, — k^ki,/k‘^)D{—k‘^). We assume that quarks 
are massless and normalize the propagator at the Euclidean point 


D{—k'^) = —d{—k^/fjj^,as)/k‘^, d{l,g‘^) = 1. 


Invariance under the RG, in the Landau gauge, leads 


d{—k‘^/lj?, as) = exp 



(-fc2) 


(5{x) 



(26) 


(27) 


The formal power series for the anomalous dimension function 7 ^, in the Landau gauge, 
is given by 


'Jv (Q^s) 


'll 

'dvr ^ IGtt^ 


7i 2 

as 


+ ...), 


(28) 


with 7 o = (13 — 4iVj/3)/2 and 71 = (531/8 — 61A^//6). 

A detailed investigation of the gluon propagator in the framework of analyticity and 
asymptotic freedom has been undertaken in Refs. 
unsubtracted spectral representation 


The gluon propagator satishes an 


D{k‘^,as,fi) 


00 



0 


pv{,^1 h) 

a — k"^ — iO 


da. 


(29) 


For a limited number of flavors < 9 in the Landau gauge the weight function satishes 


the super convergence relation / p^a^ag, p)da = 0. It was argued in Refs.|21| that For 

0 

Nf < 9 there is an renormalization invariant point K‘^{as,p) such that p^ is a negative 
measure for a > K^. We can start from formula (^). The solution for the /3-function is 
known in APT |P]. But, the anomalous dimension function 7 ^ is not determined within 
this method yet. One possible way is to assume following relation 


I v,an 





(30) 


where and are standard perturbative expressions for the renormalization group 
functions in the n-loop order, and and are corresponding quantities in APT. 
Then from formula (^), in the one-loop order, we may write 



1 /'aSlzh) V"'* = _ 1 / 1 , 1 \ 

a, ) F Tln(-F/A2) l + F/Ay 


( 31 ) 


^ (P^) is still valid in axial-type gauges, as well as, in the background field gauge. 

















with c„ = 


47r 


.^ 7 o// 3 o_ ji; jg easy to convince that expression (^) has correct analytical 
properties, since, is analytic and it does not vanish in the hnite part of the 

complex fc^-plane. Thus, function (|3lD satishes the Kallen-Lehmann representation. The 
corresponding weight function is given by 


where p+ ( p_) is positive (negative) measure, p+{cr) = Cy7rd(a) and 


(32) 


p_(a) = p_(ct) = -c. 


1 


ah? 


7n 

[i?(d]^o sin 


7o 

( ^ ^ 

— arccos 
_Po 

[R{d) [ 


In a 


+ 


\n a + a + 1 


here a = a / h? and 




In a 


\ l^ln^cr + TT^ 1 + aJ Vln^a + vr^ 


+ 


TT 


The superconvergence relation in this case reads / p-{a)da = —'ncy. For the RG invariant 

0 

scale which has been introduced in Refs.|2l[], the solution (^ID predicts trivial 

value K‘^{g,p‘^) = 0. The weight function (^) has integrable singularity at a = 0. 

For the two-loop gluon propagator similar procedure leads 




(33) 


where 

d?\a) = (a)^°/^“(l + 

it is evident that (|3^ satishes requirement of analyticity, since, the equation q;®(— fc^) = 
—dyr/do/A < 0 has no solutions in the cut complex fc^-plane. The reason is that a?2{~k‘^) 
itself satishes the Kallen-Lehmann representation with positive dehned weight function. 
One can assume the relations similar of (^) for the quark and the ghost helds too, and 
by the way, write ’’analytically improved” solutions for the corresponding propagators. 

Now we suggest diherent method for restoring analyticity and renormalization invari¬ 
ance in perturbation theory. In this case we do not accept the relation (|30D. Instead, 
we start from the perturbative expression for the gluon propagator. To the one-loop 
approximation, in the Landau gauge, we have 


d^^\—k‘^/p‘^, cxs) = 


i + StoK-S) 


'A. 

,70^5 




) - * 1>1(0 ’ 


(34) 


\ ' 4 y 11^ j 70 

here 70 is given by (|^) , ^ = h?/p^ and we have used standard perturbative 

dehnition of the A parameter h? = The ”analitization procedure” for ( [5^) leads 

(«) «s) = 4n («) «*)/4n (1) «*) (35) 
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where 


4n(w,«J = 


In M — — In f 

■wa 3 



(36) 


and u = —In order to preserve the normalization (see (PB|)) we have introduced in 
(^) the necessary nonperturbative factor. Furthermore, the invariant scale A in (^) now 
is determined according APT, since, the procedure implies the replacement Kpt Kan- 
Thus, to hrst order according APT P], we have the relation 


_ /3o _ 

dg “ otg 

47r 




1 

= - + 


where 0 = — In^. Using 


In^ 1 — ^ (p ' I — e't’ 
we hnd the anomalous dimension for the gluon held 


du 

0 l_e-(/3o/To)'/> 


02 

(l e-(/^o/T'o)'#’)2 

(70/do) 



(37) 


The expression (pSf) does not obey renormalization invariance. We can consider ([35|) as 
a known estimate and then use Eq ( 0 ) to derive the improved solution. Taking into 
account the identity 

ras{-k'^) . fln(-k^/A^) 


I3{x) 


■dx = 


7 ln(/x 2 /A 2 ) 


from Eqs. (|^) and we obtain 




70 


7^(0)d0, 


1)^ 


2„t 


Ve 


00 in(/.2/A2) t(e* - 1) (e* - 1 - t) 


dt 


(38) 


The integral in (P^D can be performed explicitly. This yields 


where 


iio/M 




^ln(-A:VA2) 


-F/A2 


I AO 


(39) 


(40) 


We see that function (^) satishes the Kallen-Lehmann analyticity for 7o/0o A 0.5. The 
corresponding havor condition is Nf < 3. For Nf > 3 the unphysical singularities appear 
in the hrst Riemann sheet. This limitation for Nf turns out to be natural. Since, the 
singularities in (|40D occurred at |fc^/A^| = 1 where the number of active quarks are just 
three. For this reason we cannot reject the solution (^Of) using arguments of analyticity 


^The author owes to D.V. Shirkov for drawing his attention to this peculiarity of the solution (p 
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Thus, we need independent criterion to select the suitable analytic solution. Recently 
Such a criterion has been suggested in Ref-ljT^, where the principle of minimality for the 
nonperturbative contributions in perturbative (ultraviolet) region has been formulated. 
Using this principle one can easily verify that the solution 


is preferable. Indeed, (|39|) 
predicts more rapid decrease of the nonperturbative contributions in ultraviolet region 
then the solution (^) provided that 66/39 < /do/bo < 2 for 0 < < 3. 


4. Conclusion 

The structure of the /3-function has been analyzed in APT. Modified expansion for 
the /3-function is written as power series in as with coefficients which depend on the QCD 
scale parameter A. The coefficients are regular at A = 0. 

We have found compact expression for the analytic two-loop effective coupling constant 
in terms of the Lambert W function. This expression is convenient to explore analytical 
properties of the solution. 

The method of APT has been extended for the gluon propagator in the Landau gauge. 
In this case, as we have seen, the nonperturbative ambiguity arises. The RG and ana- 
lyticity constraints alone are not sufficient to uniquely determine solution for the gluon 
propagator. It was shown that application of the principle of minimality for the nonper¬ 
turbative contributions in perturbative (ultraviolet) region allows us to reduce the 
ambiguity. 

The procedure does not destroy renormalizability of the theory. The only effect is 
that the renormalization constants receive hnite nonperturbative corrections. This was 
demonstrated by explicit calculations of the /3-function and the anomalous dimension 
function of the gluon held. 
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